ON THE STEADY VISCOUS FLOW OF A NONHOMOGENEOUS 

ASYMMETRIC FLUID 



FABIO VITORIANO E SILVA 

Abstract. We consider a boundary value problem for the system of equations de- 
scribing the stationary motion of a viscous nonhomogeneous asymmetric fluid in a 
bounded planar domain having a C 2 boundary. We use a stream-function formula- 
tion after the manner of N. N. Frolov [Math. Notes, 53(5-6), 650-656, 1993] in which 
the fluid density depends on the stream-function by means of another function deter- 
mined by the boundary conditions. This allows for dropping some of the equations, 
most notably the continuity equation. With a fixed point argument we show the 
existence of solutions to the resulting system. 



Introduction 

Density dependent fluids are well known and have been studied by several authors. 
Antontsev, Kazhikov and Monakhov treat in [1] an assortment of problems on density 
dependent flows of either compressible or incompressible Newtonian fluids. A more 
recent account on such problems and some improvements on results in [1] are available 
in [10] and references therein. 

The case of non-Newtonian fluids is less studied then the previous ones. Several 
constitutive laws lead to such fluids and, among them, we are particularly interested in 
the non-symmetric fluids, named micropolar fluids, introduced in [4]. Particles of these 
fluids undergo translations and rotations as well and their theory has proved to be useful 
in describing phenomena in which the structure of the fluid should be accounted for, 
e.g., blood flow in thin vessels or flows of some slurries and polymeric fluids, see [11]. 

Some authors studied evolutionary density dependent flows of micropolar fluids, such 
as the papers [2, 3, 6] whereas, to the best of our knowledge, basic results on flows in 
a stationary regime are still lacking. Our aim in this paper is, therefore, to extend the 
result in [5] placing the theory of micropolar fluids in a similar level of knowledge as the 
theory of the standard Navier-Stokes fluids. 
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Frolov addressed in [5] the stationary 2D flow of an incompressible inhomogeneous 
fluid resorting to a judicious stream-function formulation. His formulation has been 
successfully adapted to solve new problems as the mixing of two fluids having different 
(and discontinuous) densities [13] and to the boundary control of steady state flow of a 
viscous incompressible nonhomogeneous fluid [8]. 

The field equations of the model, in a steady state regime, form the following system 



~{c a + c d )Aw + pj(v ■ S7)w - (co - c a + c rf )V(V • w) + 4p r w = 2// r V x v + pg 

which we assume to hold in a bounded planar domain f2, having a C 2 boundary <9f2, 
subject to the following boundary conditions 



with r C dQ, being a connected arc on which vq ■ n < 0, that is, T is the part of d£l 
where the fluid flows inward. 

The equations (1) represent conservation of linear momentum, the continuity equa- 
tion, the incompressibility and the conservation of angular momentum respectively. The 
unknowns are p, the denstity; v and w, the fields of velocity and rotation of particles 
and the pressure, p. The fields / and g are, respectively, given external sources of 
linear and angular momenta densities whereas p, p r , Co, c a , Cd, J are positive constants 
characterizing the medium and also satisfying Co > c a + c d - 

Notice that this model contains the incompressible, density dependent Navier-Stokes 
system as a particular case (p r = 0,w = 0), and it is named the nonhomogeneous 
micropolar fluid model. More details on the derivation of the model as well as the 
physical meaning of the several parameters above may be found in [11, 4]. 

Below we devised a way to solve problem (l)-(2) by combining arguments used in [11, 
5] . This is a new result in theory and, although its proof follows fairly well-grooved lines, 
we believe it may serve as a first step towards a better understanding of similar flows 
in different geometries such as infinite or semi-infinite pipes. 

Begining with Frolov's stream-function approach allows us to drop the continuity 
equation and the boundary condition on p, cf. Section 2. Then we consider an auxiliary 
linear problem, for the rotational velocity w, which we solve with the Lax-Milgram 
lemma, see Section 3.1. Next we solve a problem for the translational velocity, v, by 
defining a suitable operator and using the Leray-Schauder principle, cf. Section 3.2. 
Thus we get a pair v, w as a weak solution to our problem. We conclude our paper with 
some remarks concerning the solution so obtained. 



(1) 



-(p + p r )Av + p(v ■ V)v + Vp = 2p r V xw + pf 
V ■ (pv) = 0, V • v = 0, 



(2) 



p = po > on T, 
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1. Notations and statement of the main result 

We shall now introduce some notations and clarify what is meant by (1) to hold in a 
planar domain. Given u= (ui,u 2 ,u 3 ),v = (vi,v 2 ,v 3 ), we denote [(u-V)v]j = u-Vvj = 
J2k u k^x k Vj, j = 1,2,3. Loosely speaking, a planar flow may be sought as a "slice" of 
a 3D one, that is, the flow takes place in a cross-section 23 = const, of a 3D domain. 
Thus, the functions involved are assumed to be independent of the 23 variable and the 
axes of rotation of the particles of the fluid are assumed to be perpendicular to the 
plane of the flow. This way, we regard 

v = {v 1 (x 1 ,x 2 ),v 2 (x 1 ,x 2 ),0),p = p(xi,x 2 ), w = (0,0,^3(21,22)), 

/ = {fi(xi,x 2 ),f 2 (x 1 ,x 2 ),0),g = (0,0,53(21,22)), 

and write V -1 ^ = (—d X2 4>, d Xl ip) and V x (</>i, <fe) = d Xl c/) 2 — d X2 (f)i so that the system (1) 
may be written componentwise as 

+ p r )Avj + pv ■ Vvj + d Xj p = (-l) j ~ 1 2fj, r d Xj w 3 + pfj, j = 1,2 

(3) V • (pv) = 0, V • v = 

~{c a + c d )Aw 3 + pjv ■ Vw 3 + 4p r w 3 = 2p r V xv + pg 3l in tt. 

From now on, we shall adopt j = 1, a = p + p r and K = c a + to shorten the equations. 

We use standard notations regarding Sobolev spaces modelled in L q (Q), W k ' q {Q), k > 
0, q > 1, and their norms || • ||^fc,<j • The same goes to the trace spaces, W h ~ 1 / q ' q (dQ), k > 
0, q > 1, and their norms || • \\ W k-i/ q , q - For q = 2, we write W k ' 2 (tt) = H k (Q) and 
W l l 2 > 2 {dVL) = HV 2 (dn), k > 0, as usual. 

By V we denote the set of divergence- free vector fields 92 = (ipi,<p 2 ) such that ip € 
C^°(fl), V is the closure of V in the i? 1 -norm and H = {<p £ H 1 \ V • p = in Q}. 
By C m, P(£l) we denote the set of all m times continuously differentiable functions in ft 
whose m-th order derivatives are Holder continuous with exponent f3 € (0, 1). As above 
|| • stands for the C ' 13 (Q)-norm.. 

Our main result then reads 

Theorem 1. Let f,g € L 2 (Q), p € C°>P(T), < /3 < 1, and v ,w € H 1 / 2 ^) 
be given satisfying (2). There exists a weak solution p G C°' a (Vl),a < (3,v £ H, 
w £ H 1 ^), of system (3) in the sense of Definition 1 below, provided that p,K are 
sufficiently large so that min{/i,2K} > C||?7||loo 1 1 -tt?o 1 1 ^1 /a - 

2. Weak formulation of the problem 

For a given a divergence- free vector field v = (v\, v 2 ) in Q there exists <f> : ft — > R such 
that v = V^cf). In addition, denoting by r, n the unit tangent and outward normal fields 
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on dQ, and bearing in mind conditions (2), the assumption V <j) = on T amounts to 

Thus, the boundary values of 4> may be obtained upon integration, with respect to 
the arc length, from a point x G T, 4>{x) = — ff Vo • nds, x G T. Moreover <fi G 
H^ 2 (T) C C(r), as v ■ n G H 1 ' 2 ^). From v ■ n < on r we see that d> is increasing 
on r and we may speak of its inverse cf)^ 1 : <fi(T) C M. — > T. We may then define 
v(y) = Po( ( P~ 1 (y))j y ^ <fi(T) C M and extend it to K as a strictly positive function 
rj G C 0,/3 (IR),/3 < 1, such that ij(tp(x)) = po(x),x G T, whenever V^V = on T. 

We fix this rf meeting the above requirements. For sufficiently smooth 77, ij) we have 
V • [ry(V')V- L V] = ■ V^-if) = 0, in 0. A weak version of it may be proved by 

regularizing r] if its assumed to be merely continuous. From all the above facts, we 
assume that the density has the form p = r)(ip), where V^V = v, whence, the continuity 
equation (3)2 may be dropped as well as the boundary condition (2)i on p, see [5]. 

We follow [8] and denote by iV : H — >■ H 2 (Q) the continuous operator assigning to 
each divergence-free vector field, u, in f2 its stream-function ip = Nu. Actually such a 
stream-function is determined up to an arbitrary additive constant which we take to be 
zero with no loss of generality, see e.g. [14, Lemma 2.5, Chapter 1] or [7, Theorem 4]. 

Let vq,wo G H l / 2 (dQ) satisfy (2)3. From the trace theorem follows the existence of 
b G H 1 ^) with 

(4) 6|qq = wq and the estimate \\b\\ H i < C||i/;o||#i/2, 

holds for some absolute constant C. Also, for each 5 > fixed, there exists the so-called 
Leray-Hopf extension of vq, that is, a vector field a G H 2 (VL) satisfying, 

S\qq = vq, / a 2 ■ ip 2 dx < 5 2 |V(/9| 2 dx, ip G V 

(5) Jn Jn 

V ■ a = 0, in ft, a(x) = 0, d(x, dfl) > e, for a fixed e > 0. 

For a construction of the above cut-off functions, the reader is referred to [9]. From this 
point on we shall omit the domain of integration since no boundary integrals appear in 
the subsequent computations. 
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Next we introduce new unknowns u = v — a, w = w — b, u G V and w G Hq(Q), 
satisfying the following system of equations in ft 

-aAu + p[(u ■ V)u + (a ■ V)u + (u • V)a] + Vp = -2p r V ± w + pf + F 

V • (pu) = 0, V ■ u = 0, 

( 6 ) 

—kAw + p[u ■ Vu> + u ■ Vb + a • Vw\ + 4p r w = 2p r V x u + pg + G, 

u,w = 0, on c?Sl, 

with F = F(p, a, b) := aAa - p(a • V)a — 2p r V ± b and G = G(p, a, b) := nAb - p(a ■ 
V)6 — Ap r b + 2/i r V x a. 

To solve the boundary value problem (6) we look for p = i](N(u + a)) and set 
F = F(u,a,b) := F(r](N[u + a]),a,b) and G = G(u,a,b) := G(rj(N[u + a],a,b). From 
the imbedding H 2 (Q) C C°' e (Q),9 G (0,1), and r? G C^R) we see p G C°> a {Tl), for 
a = £0 < p. 

Definition 1. We call a pair u G V, w G -ffo(^) « weaA: solution of system (6) i/ i/ie 
integral identities 

a j Vu-X7<pdx = J rj(N(u + a))[u ■ Vtp ■ u + a ■ V<p ■ u + u ■ V99 • a] dx 

— 2p r / V^w ■ ip dx + / n(N(u + a))f ■ (p dx + / Fipdx 

(7) , 7 " 7 ' 7 

k / VwV^dx= j r)(N(u + a))[(u- V£)w + (u- V£)6 + (a- V£)w]dx 

- 2pi r J (2w - V x u)£dx + y ?](N(u + a))g£dx + J G£dx 

hold for all (p G V, £ G -ffg (^) anc ^ ^ ^ as above. 

It is readily seen that a weak solution u, w of system (6) in the sense of the above 
definition yields a weak solution v = u + a, w = w + b to the original problem (l)-(2). 
Indeed, recalling (4) and (5), we see that the boundary conditions (2)2 hold in the sense 
of traces. Moreover the integral identities (7) for u = v — a, w = w — b imply analogous 
ones for v, w. 

It is worth remarking that the recovery of the pressure is a standard matter and it 
follows from the knowledge of p,v,w, see e.g. [14]. To be precise, we state it below as 

Theorem 2. Let p G C 0>/3 (n),0 < /3 < l,v € H,w €. H 1 ^) be a weak solution of 
the problem (l)-(2). Then there is a p G L 2 (fl), such that Vp G L 2 (fl) and for all 
ip G C^(n) it holds that 

a J Vv ■ Vip dx — j p((v ■ V)"0 • v — f • if) dx + 2p r J \/ ± w ■ ipdx = j pV ■ ip dx. 



6 fabio vitoriano e silva 

3. Proof of Theorem 1 
We split this proof in three steps. 



3.1. Auxiliary problem. In this section we consider the following auxiliary problem: 
(A) given »6H, find u> G H$ (Q) such that the identity 

-kAw + r](Nv)v ■ Vw + 4fi r w = 2fi r X7 x v + i](Nv)g + G(v, b), in ft 

holds in the sense of distributions, with G(v, b) = nAb — r](Nv)v ■ V6 — \[i r b. Existence 
of an unique w G Hq(Q.) solving problem (A) follows from the Lax-Milgram lemma. In 
fact, the problem (A) amounts to find w G Hq(Q) such that 

(8) B v \w,£] = J (2n r V x v + rj{Nv)g + G(v,b)^dx, for all f G 

and B v : Hq(Q) x Hq(Q) ->■ E being defined as 

= y [kVx • Ve + V (Nv)v ■ V X £ + 4/vx£] <*c, for all X , £ € ^ (O) . 

According to standard estimates B v is readily seen to be continuous. 

It is also coercive as for all v G H, w G Hq(Q) we have J i](Nv)v ■ Vwwdx = 0, 
whence 

B v \w,w] = k\\ Vu>||^2 + 4/x r ||w;||^2 > min{K, 4// r }||u/||^i, 

for all w G Hq(Q). 

In addition, the right-hand-side of (8) is a continuous form in Hq(Q), 



J (2fi r V x v + r](Nv)g + G(v, b)^dx 



<C\\Z\\hi, 



for all £ G Hq(Q), with a constant C depending on c a ,Cd,f^ r , \\v\\l°°, IIs'IIl 2 ) INo||#i/2, 
\\v\\jji and 0. The Lax-Milgram lemma assures the existence of an unique w G Hq(Q) 
solving problem (A). 

We point out for future reference that the following estimate holds for w: 

«||Vu?||^2 + 4/i r ||7£||| 2 < 2/x r ||V7;||£,2||u;|| £ 2 

(9) + ||??||l°°||5'IIl 2 II^IIl 2 + CU^IlL^lluJollHi/all^ll^lNIlL 2 

+ CK||tyo||// 1 /2||Vw|| L 2 + 4:fi r C\\w \\ H l/2 \\w\\ L 2. 

3.2. Problem for u. Our goal is to obtain u G V solving (11) as a fixed point of 
the operator A, to be defined below. To this end we shall apply the Leray-Schauder 
principle, which requires the operator A to be completely continuous and also that every 
possible solutions of u = XAu, A G [0, 1], are uniformly bounded, see [11, 9]. 



STEADY FLOW OF NONHOMOGENEOUS ASYMMETRIC FLUID 



Define A : V — > V as follows: for u G V, let u; G #o(^) denote the solution of 
problem (A) corresponding to v = u + a and consider Au given by the identity 



(10) 



a j V(Au) -V(pdx = J rj(N[u + a])[(u ■ V)<p -u + (u-V)<p ■ a + (a ■ V)ip ■ u] dx 

- J {^HrV^w- rj(N[u + a])/) -<pdx + J F-<pdx, 

for all (p G V. The well-definiteness of *4 follows from Riesz theorem. Indeed, the rigth- 
hand side of (10) is a continuous form on V owing to the properties of a, cf. (5), and 
some elementary estimates. 

Now we check A is completely continuous. For G V, i = 1, 2 let us define = 
{f + a G H and consider w l G iTg (^) the solutions of problem (A) corresponding to v 1 . 
Then for ip G V, i = 1,2, 

a y V(«4if) • V<^cte = ^ r/(A^[if + a])(^ • V)v? • if 

+ y ^(JV^ + a|)[(« i • V)<^-a + (a- V)y» ■ «*] ds 

- y ^2/i r V ± u; i - r}{N[i? + a])/) ■ ipdx + j Ffipdx, 

with Ft = F(ri(N(v* + a), a, b). 

Subtracting these two identities we get 



(11) 
for 



a y V{Au 2 - Au 1 ) -Vipdx = h+ I 2 + h + h + h + h, 

h = J (v(N[u 2 + a]){u 2 ■ V)(p ■ u 2 - 7]{N[u l + a]){u l ■ V)(p ■ u 1 ) dx 

1 2 = y (r](N[u 2 + a])(u 2 • V)v? • a - rj{N[u l + a])^ 1 • V)<^ • aj dx 

1 3 = J (^(N[u 2 + a])(a ■ V)tp ■ u 2 - rj{N[u l + a\)(a ■ V)<p • m 1 ) dx 

14 = — y 2/j ir \7 ± (w 2 — w 1 ) ■ (pdx 

= + y 



/ 5 
/ 6 



+ d])-??(ivp 1 + d]))]/-^c?x 



= ~y r/(AT[n 2 + a])) - ??(iV[u 1 + a])) (a-V)a-(pdx. 
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We now show it is possible to bound each of the Ik, k = 1, . . . , 6, by a constant times 
||u 2 — - u 1 || i 4, proceeding as follows. First notice that 

h = J [v(N[u 2 + a}) - ^(Nlu 1 + a])) (u 2 ■ V)<^ • u 2 dx 

+ J ??(iV[u 1 + a]){{u 2 ■ V)ip ■ u 2 - (u 1 • V)<p ■ u 1 ) dx 

= J (v(N[u 2 + a]) - //(iVp 1 + a])) (u 2 ■ V)tp ■ u 2 dx 

+ JriiNlu 1 + a])|[(n 2 - u 1 ) ■ V]<p ■ u 2 + (u 1 • V)ip ■ (u 2 - u 1 )} dx = I u + Ji 2 . 
Holder, Young and imbedding inequalities imply 

\hi\ < Mo j \N{u 2 -u^r^uld^u^dx 



j,k j,k 

^Cll^llc-llV^ll^lliV^-n^ll^ll^llis 
= C||7 ? ||^||V^|| L2 N 2 -n 1 ||24«N 2 |||8, 

I/12I < ll^ii^iiv^il^di^ii^ + iin 1 !!^)!^ 2 -^ 1 !!^. 

Hence 

j I / /~>ll„ll IIV7.-II IL-?2 „-?l| 

L 4a II u IIl 8 



|/l| < C\\T1\\C° II V^|| L 2 ||« 2 -U 1 ||24 a ||« 2 " 2 



(12) 

+ IMM|Vv9|| L2 (||£ 2 || L 4 + ||^|| L 4)||^ 2 - n 1 !!^. 

Likewise we get the bounds 

\h\ < ll^llL-llV^II^H^-^ll^llall^ 

+ Cll^llc^ll^llLsllallLsllV^II^Hu 2 - n 1 ||24 a 

|/ 3 | < IMiLHIv^lMl^-^lbiiaib 

+ C|^llc«|l^llL8||a|| L 8||V^|| L2 ||u 2 - u 1 ^ 
\h\<\\v\\c^-u^\\f\\L4<Ph* 

|/ 6 | ^cii^ii^Plllsllv^n^iin 2 -^ 1 !!^. 

The term J4 requires estimating \\w 2 — w}\\ H i in terms of \\u 2 — u 1 ^, a task we now 
perform. As w l , i = 1,2 solve problem (A) we find that the following identity 

J kV(w 2 - w 1 ) ■ Vip + 4n r (w 2 - w})tp dx = Ji + J 2 + J 3 , 
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holds for all tp € Hq(Q), where 

Ji = J (r](N[u 2 + a])(u 2 + a) ■ Vw 2 - r](N[u l + a])(u 1 + a) ■ Vw^ipdx 

J 2 = J (2fi r V x (u 2 - u 1 ) + [r](N[u 2 + a])(u 2 + a) - ? ? (iV[u 1 + a]){v} + a)]g^)ipdx 

Ja = J (r]{N[u 2 + a])u 2 - r,{N[u l + a])^ 1 ) • Vb^dx. 

By arguing as in the estimations (12)-(13) we find 
\Ji\ < c(jMlc«N 2 - ^ll^^ll^ 1 + a\\ L & + II^HzHI^ 2 - ^IIl 4 ) IIV^IIlzII^IU 1 
|J 2 | < C(2^ r meas(0) 1/4 ||n 2 - u l \\ L ± + |M|o \\g\\ L 2 \\u 2 - u^Jj \\^\\ H i 

IJ3I < ^(ll^llcHI^ 2 - ^ll^^ll^ 1 + o|| L 8 + ||??||l-»||m 2 - ^Wl^j \\vb\\ L 2\\ip\\ H i. 

Taking ip = w 2 — w 1 and invoking the boundedness of we get 
\\w 2 - w^Whi < c(\\Vw 2 \\ L 2 + ||V6|| L 2^|||r7||c«||'u 1 + a|| L 8||n 2 - u 1 



\J1 „-?l|ia 

II 4 



+ \\u 2 - w 1 1 [ ^4 j> + (2^ r meas(0) 1 / 4 + \\rj\\c<* \\9\\l 2 )\\u 2 - ^Hl 4 , 

for some constant C > depending on K,[i r ,Q. At last, bearing in mind the triangle 
inequality, we find 

\h\ < 2/i r c(||V^ 2 || L2 + ||V6|| L2 )||r ? || c «(||u 1 || L 8 + ||a|| L 8)||n 2 - uX^WfWm 
(14) v ' 

+ 2^ r C(2 / u r meas(0) 1/4 + \r}\ c <* IIs'IIl 2 )!!^ 2 ~ ^IMMItf 1 - 

Collecting inequalities (12)-(14) for ip = Au 2 — Au 1 , together with \\u 2 — u 1 \\l4 < 1 
and in view of (4) and (5), we conclude 

a\\V(Au 2 - Au^Whi < \\u 2 -u l \\ a LA x {ciHIc-o ||n 2 || L 8(||{i 2 || L 8 + C\\vo\\ H i /2 ) 

+ CllVoll^dl^HiS + C\\V \\ H l/2) + ||/1| L 2 + || 5 || L2 

+ (ll^hs + C\\v \\ H i/2)(\\Vw 2 \\ L 2 + C\\w \\ H i/2) 

+ \\v\\l°° II^ 2 IIl 4 + II^IIl 4 + 2c 1 1 v \ \ h 1/2 

+ (ll^hs + C\\v \\ H i/2)(\\Vw 2 \\ L 2 + C\\w \\ Hl/2 ) 
+ 2/i r meas(^) 1 / 4 ] } . 

Thus, from the compactness of the imbedding H 1 ^ L 4 and the inequality above, a 
weakly convergent sequence in V is mapped by A into a strongly convergent sequence in 
L 4 . It remains to show that all possible solutions of the equation u = \Au, A G [0,1], are 
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uniformly bounded. As A = implies u = 0, we suppose A > and replace Au = u/A, 
ip = u in the equation (10). Performing estimations similar to the previous we find 
(15) 

u-adx +2/i r (|H| i2 + ||6||z,2)||Vtt|| L 2 + ||?7||l°° ||/||l 2 NIl 2 - 



& ,,„^,,9 



Notice that the term 2// r ||u;||£2 1| Vt?||x,2 itself is another quadratic term in u which we 
handle as follows. In view of (9), with v = u + a, we find 
(16) _ 

^||Vu7||| 2 +4// r ||«;||| ( 2 < 2/i r ||V(tT + a)|| L 2||«;||x,2 + \\r]\\ L <x> \\g\\ L 2 \\w\\ L 2 

+ C||?7||Loo||-u; || Hl/2 ||n + a|| L 4||u;|| L 4 + C/c||«;o|| H i/2 1| Vw|| L 2 
+ 4/i r ||-u;o|| H i/ 2 ||u)|| L 2. 

Aided by Young and Holder inequalities we get, by summing up (15) and (16), 



A 



\\Vu\\l, + 4Vw\\t 2 <\\ V \\ L ~ 



J (u ■ V)u ■ adx 

+ fi r \\Vu\\ 2 L 2 + C||t7||loo ||«||x,4||it7o|| jff i/2||u;|| L 4 

+ 2/x r ||6|| i2 ||Vu|| i 2 + ||^||l-»||/||l2||u|| L 2 

+ 2^ r ||Va|| L 2||u;|| L 2 + \\r)\\L°° WoWl 2 \\"ML\\l 2 

+ Cll^lk 00 1 1 « 1 1 l 4 I No 1 1 #1/2 \\w\\l 4 + C k I|\7w||l2||u>o||#i/2 

+ 4:H r \\w \\ H i/2\\w\\ L 2. 

Next, requiring S > in (5), to be such that 5||r/||x,oo < fi/2 and estimating 

C / \ 

II^IU oo ll«lll,4||«'0||Hi/2||«i|| £ 4 < — H^Hioo ||w Hjyi/a (J|Vu||| 2 + HVwIllaJ 
we ultimately get, using (4) and (5), 

c c 

(.2X ~ ^\\L°°\\wo\\m/z)\Nu\\ 2 L2 + ( K ~ 2-|I j 7IU °II^o|Ihi/2)||V2«|| l2 

< 2fl r C\\w \\ H l/2\\Vu\\ L 2 + |M|l°°|I/IIl 2 NIl 2 + 2/i r C||t/ || H l/2||lt)|| L 2 

+ IM|l°°||5Hl 2 IMU 2 +C||77||l«.||t;o|| h i/2||«;o||hV2||^||l4 
+ C«||V«;||£2||t(;o|| H i/2 + 4^ r C||t(;o||/fi/2||w||i,2- 
Therefore, demanding [i and k to be large enough so that 

min{^,, 2k) > C\\t]\\l°° \\wo\\ h i/2, 

we may conclude the uniform boundedness on the norms of all possible solutions of 
u = XAu, A G [0,1], with respect to the parameter A. From this and previous steps 
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we conclude that the Leray-Schauder principle applies and that problem (6) has a weak 
solution. □ 

4. Concluding remarks 

We have tacitly assumed dfl to consist of a single component, that is, ft C M 2 to be 
a simply connected open set. A clue on how the arguments should be modified to cope 
with a non-simp ly connected ft may be found in [5]. 

The presence of the quadratic term brought up by the equation (1)4, for the rotational 
field w, required us to demand the viscosities fi, k to be sufficiently large compared to 
data and this somewhat contrasts with previous results by Frolov [5] and Santos [13]. 

As shown in [5, cf. Theorems 2 and 3], also [11, Chapter 2], the regularity of p,v,w,p 
may be improved by increasing those of dfl, po,vq,wq, f and g. An usual bootstrap 
argument should combine regularity results for the Stokes problem and Necas' results 
on strongly elliptic systems of second order [12]. 

We also notice that the uniqueness of the above solution deserves further investigation. 
Indeed we could not benefit from neither the steady-state continuity equation (1)2 nor 
the particular form of the density to derive the required estimates. 
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